Notes on coherent backscattering from a random 

potential 



Herbert SpohE0 

Zentrum Mathematik and Physik Department, TU Miinchen, 
D - 85747 Garching, Boltzmannstr. 3, Germany 



Abstract. We consider the quantum scattering from a random potential of strength 
A 1 / 2 and with a support on the scale of the mean free path, which is of order 
A -1 . On the basis of maximally crossed diagrams we provide a concise formula for 
the backscattering rate in terms of the Green's function for the kinetic Boltzmann 
equation. We briefly discuss the extension to wave scattering. 
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1 Introduction 



The motion of a single quantum particle in a weak random potential of strength 
A 1 / 2 is well approximated by a kinetic equation. More precisely, in the limit of small 
A, the disorder averaged Wigner function W is governed by the transport equation 

^W(q,p,t)) = -V p cu(p)-V q W(q,p,t) (1.1) 

+p(q) [ rfp , 27r5(^(p)-^(y))(27r)- 3 / 2 %-p')(^(g,p',t)-^(g,p,t)). 



Here w(p) is the kinetic energy (= dispersion relation), u>(p) = p 2 /2m for a non- 
relativistic particle. 9 is the Fourier transform of the translation invariant part of 
the covariance of the random potential and p(q) 1 ^ 2 is the spatially varying strength 
of the potential. The mean free time and the mean free path are both of order 
A -1 and Eq. (jl.ip is written on that scale. As a side remark, the Wigner function 
integrated against a smooth test function is self-averaging [T] . Thus (11.11) holds even 
for a typical realization of the random potential. 

(11. ip can be thought of as arising from the stochastic motion of a fictitious 
classical particle: the particle moves along a straight line according to its velocity 
and randomly changes its momentum, respecting energy conservation, with a rate 
given through the covariance of the potential. In this sense the wave nature of the 
quantum motion is completely lost as A — > 0. However, as discovered by Langer and 
Neal p], even for small A the wave character is still visible in the fine structure of 
the average scattering rate. We imagine an incoming plane wave with wave number 
k scattering off the random potential and denote by 0b(&, k') the scattering rate 
from the incoming wave number k to the outgoing wave number k! according to the 
transport equation (II. ip and by (<J\(k, k'))y the scattering rate for the Schrodinger 
equation averaged over the random potential. Then 

lim X 2 (a x {k, k')) v = a B {k,k') loxk'^±k. (1.2) 

A — >0 

Note that the total cross section of the scattering potential is of order A~ 2 . Therefore 
in (11.21) one needs to balance by the prefactor A 2 . For exact backscattering one finds 
that 

lim X 2 (a x (k, -k)) v = 2a B (k, -k) , (1.3) 

A — »0 

where the tilde indicates that from <tb the paths with only a single scattering event 
have to be omitted. In fact, the enhanced backscattering consists of a narrow peak 
of width A centered at k! = —k. As to be discussed, one finds 

lim \ 2 (a\(k, -k + \k)) v = 0b ac k,fc(«) . (1.4) 

A — *0 

Properties (11.21) and (11.41) can be seen from an expansion of the scattering ampli- 
tude in powers of the potential V. Summing all ladder diagrams yields (II. 2p . while 
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(11.41) results from the maximally crossed diagrams. Of course, to actually prove the 
limits (II. 2p . (11.41) one would have to establish that the contribution of the very many 
remaining diagrams vanishes as A — > 0. This task will not be addressed in my notes. 
But I point out to the reader that very complete and precise estimates have been 
achieved in the recent work of Erdos, Salmhofer, and Yau [3l Hj. These notes are 
written, at least in part, with the aim to encourage a similar kind of analysis for the 
scattering rate. 

Coherent backscattering according to (II. 4p has been confirmed experimentally 
including the shape of the peak [51 [6] . The experiment is done for light scattering, 
since light is so much easier to manipulate than electrons. In ( 11.41) averaging over 
disorder is required. In a single scan no peak can be disentangled. Thus one either 
repeats the experiment many times using samples of disordered glass or, more ele- 
gantly, scatters from a turbid solution. Then the Brownian motion of the suspended 
particles provides the averaging for free. Excellent theoretical texts are available 
[TJIEIE]- They are striving for even more refined information, as e.g. the variance of 
the conductance fluctuations and the statistics of speckle patterns [7J. Does there 
remain then anything to be done on the level of (11.41) ? 

In fact, I was looking for a concise expression of <7back,fc in terms of the transition 
probability resulting from Eq. (II. ip and could not find it in the literature. My result 
is given at the beginning of Section The required computation I find sufficiently 
illuminating to be put into written form. 

In the literature (as far as I checked) the conventional approach is to start from 
the average of the square of the Green's function and then to extract from it the 
scattering. We proceed in a way which looks more systematic to us. For every 
realization of the disordered medium the scattering rate is given through stationary 
scattering theory, at least in principle. The average scattering rate is then expanded 
in the disorder strength. 

The standard Born expansion for the scattering amplitude is covered in Section 
|2] including the average over the random potential. The scattering theory for (ll.ip is 
explained in Section [31 The most lengthy part of our contribution is the summation 
of the ladder diagrams, while the maximally crossed diagrams then easily follow, see 
Sections H] and [5j We add the modifications required when the Schrodinger equation 
is replaced by a wave equation and close with the diffusion approximation to (II. 4p . 

2 Average scattering rate for electrons 

Because it is slightly simpler and more familiar, we first discuss the stationary 
scattering theory for an electron moving in a random potential. The physically 
more relevant wave scattering will require only minor modifications, see Section 
We consider the hamiltonian 

H = ~^A + A 1/2 \/(x), A>0, (2.1) 
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acting on the Hilbert space 7i = L 2 (M. 3 ,dx). In Fourier space the Schrodinger 
equation then reads 

i^(k, t) = u(k)$(k, t) + (2tt)- 3/2 A 1/2 J dk l V{k - k x )$(k u t) , (2.2) 

where we introduced the dispersion relation u(k), 

u(k) = \k 2 (2.3) 

for a nonrelativistic particle of mass 1. The formalism is such that general disper- 
sion relations are allowed. V(x) is a Gaussian random potential with zero mean 
(V(x))v = 0, (-)v denoting the average over V. The statistics of V(x) is translation 
invariant locally but V(x) vanishes outside some bounded region. To achieve it we 
introduce the smooth shape function p : R 3 — > R, p > 0, such that the support of p, 
supp p, is a bounded set. Then the covariance for V(x) is given by 

(V(x)V(y)) v = p{\x) l ' 2 d{x - y) P {\y) 1/2 . (2.4) 

6(x) = 9{—x) and 9 is assumed to have a rapid decay. Its Fourier transform satisfies 
6(k) > 0, so to have a positive definite covariance. 

In Fourier space the covariance (12.41) is peaked near k' = —k and we write 

(V(k)V(—k + \kx))v 

= (2tt)- 3 / 2 I dg^g + ^X-^iX^gyX-^X^g + h) 



9(k)X- 3 (2n)- 3 / 2 J dgpV%g)*pV 2 (g + h) 

9(k)X- 3 (2n)-V 2 [ dxp(x)e- ikl - x = ^X^pih) . (2.5) 



Therefore, in approximation, we set 

(V(k)V(k')) v = 9(k)px(k + k') , (2.6) 

where f)\{k) = A _3 p(A _1 A;). Below we will always use (12.61) . with symmetry and 
positivity restored through the 0(A) corrections. 

We consider a weak potential, A C 1. Then the mean free path of the electron 
is C(A _1 ). The shape function p enforces V to vanish outside a bounded region on 
the scale A -1 . Thus the total cross section is of order A~ 2 . 

For almost every realization of V the hamiltonian H has a well-defined unitary 
scattering matrix S\, see [10] Section XL 6 for a discussion of stationary scattering 
theory. In momentum space S\ has the kernel 

S x (k, k') = 5{k - k') - 2mT x {k, k')5{u{k) - u{k')) . (2.7) 
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Let us set (pk(%) = 3 / 2 e* fc ' z . Then the T-matrix is defined through 
T\(k, k') = lim> fc , [\ 1/2 V - \ 1/2 V(H - (E + ie^X^V)^ , 

e— >0 

E = uj{k) = u(k') , (2.8) 

with (•, denoting the scalar product of TC. The scattering rate a\{k, k') from k 
to k! is the square of the T-matrix. Since in (12.71) k refers to the far future and k' 
to the far past, one obtains 

<j\(k, k') = (27r) 3 \T x {k', k)\ 2 2ix5{uj{k) - uj(k')) (2.9) 

for k ^ k'. Recall that a\(k, k') is random through the dependence on V. 

In these notes we deal only with the average scattering rate, (cr\(k, k'))y, for 
small A. The variance of a\ is also of great interest and has been studied both 
through experiments and theoretically [7]. 

To be able to perform the Gaussian average we expand the resolvent of ( 12.81) into 
the Born series as 

oo 

T\(k', k) = ^(-l)-(^, \ l l 2 V{G E+ \ l ' 2 Vr Vk ) n . (2.10) 

n=0 

Here Ge+ is the free Green's function, 

G E+ = lim (H — (E + te))' 1 , H = - . (2.11) 
Inserting ( 12.101) in (12.91) one obtains 

oo 

(a x (k, k')) v = (2tt) 3 (| Y^{-l) n (<p ¥i A 1 /V(G B+ A 1 /V)> fc ) w | 2 ) v 27r5(a;(A ; ) - u(k')) . 

n=0 

(2.12) 

At this stage the average over V can be carried out explicitly. The result is most 
concisely expressed diagrammatically. 

k] k-> A;; 



(a) (b) 



k'- 



91 (c) (d) 

Figure 1: A diagram and its basic building blocks. 

We draw two horizontal lines, see Fig. 1, the top one is directed to the right 
and the bottom one towards the left. They have external momenta k', k. The 
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top line carries n+ vertices, n + = 1,2,..., and corresponding internal momenta 
ki, . . . , k n i- Each bond, see (a), with momentum k carries the Green's function 
GE+{k) = (w(k) — (E + is))' 1 . Each vertex, see (b), with momenta k±, k 2 car- 
ries the potential X 1 ^ 2 (2n)~ 3 ^ 2 V(ki — k 2 ). Correspondingly, the bottom line carries 
n_ vertices, n_ = 1,2,..., and internal momentum gi, . . . , g n _-\. Each bond with 
momenta g carries the Green's function GE_(g) = GE + {g)* and each vertex with 
momenta gi, g 2 carries the potential X 1 ^ 2 (2tt)^ 3 ^ 2 V (g\ — g 2 ) ■ The resulting expression 
is integrated over all internal momenta k±, . . . , fc n+ -i, <?i, • • • , <7n_-i- The Gaussian 
average generates a sum over all pairings, provided n + + n_ is even. Otherwise the 
average vanishes. A pairing is indicated by a wavy line, see (c). 

For obvious reasons the building block (d) is called a gate. The sum over all 
gates can be performed thereby modifying the free Green's function (or propagator) 
to the effective medium Green's function as 

oo 

n-^-r— ■ 

Thus the diagrams remain as in Fig. 1 only the light lines are replaced by thick lines. 
Because of the spatial cut-off in the potential the effective medium Green's function 
is not as explicit as in the translation invariant case. 

In the following we study two particular classes of diagrams, namely the ladder 
diagrams and the maximally crossed diagrams, schematically represented in Fig. 2. 




Figure 2: A ladder diagram and a maximally crossed diagram with effective medium 
propagators. 



We denote by Ii a d{k', k; A) the sum over all ladder diagrams with external momenta 
k' and k and by I max (k', k; A) the sum over all maximally crossed diagrams with 
external momenta k' and k. 



3 Scattering rate from the Boltzmann equation 

For a weak random potential one can approximate the true wave dynamics by the 
stochastic motion of a fictitious classical particle. It has position q e M 3 , momentum 
p 6 I 3 , and kinetic energy uj{p). It moves with constant momentum and changes 
its momentum at random times subject to the constraint of constant energy. More 
precisely, given the current location q and momentum p, in the short time interval 
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dt the probability for the momentum to be scattered into the volume element dp' is 
given by 

p(q)(2n)- 3 / 2 6(p - p')27T5(u(p) - u{p')) (3.1) 

independently for each small time interval. The corresponding stochastic process is 
determined by the backward generator L, which is defined through 

Lf(q,p) = V p uj(p)-V q f(q,p) (3.2) 
+p(q) [ rfp'27r5(^(p)-^(j 5 / ))(27r)- 3 / 2 %-p')(/(g,p , )-/(g,P)) 



as a linear operator acting on functions / : M 3 x IR 3 — > C. The transition probability 
from (q,p) to dq'dp' in time t, t > 0, is then given by the integral kernel of e Lt , 
denoted by 

e Lt (q,p\dq'dp') , £>0. (3.3) 

Clearly, outside of suppp the collision operator vanishes and the particle moves freely. 

The stationary scattering situation is easily modelled. We choose w G M. 3 away 
from suppp and a plane A m j through w and orthogonal to the incoming unit wave 

vector k = k/\k\ and not intersecting suppp. The initial position of the particle 
is uniformly distributed on A m j and the initial velocity equals Vu;(fc). The corre- 
sponding initial measure p^ imposes a uniform flux at momentum k through the 
plane A ?, i.e. 

Hk{dqdp) = 5(p — k)5((x — w) ■ k)\k ■ Vu(p)\dqdp . (3.4) 
The Boltzmann scattering rate from k to k' is then given by 

2n5(u(k) -u{k'))a B {k,k') = hm J J p k (dqdp)e Lt (q,p\dq'dp')S(p' - k') (3.5) 

for k ^ k! . For k = k' there is a 5-contribution coming from those paths which 
pass the potential without any scattering. In our set-up this (^-function has infinite 
weight. More meaningful would be to restrict fik{dqdp) to those initial conditions 
for which the particle actually hits suppp. But the strict forward scattering k' = k 
is of no concern to us here, anyhow. The right side of (13. 5p is proportional to 
5(u(k) — u)(k')). It is convenient to remove this overall factor. On the energy shell, 
o"b(&? k') is uniquely defined and depends smoothly on its arguments. 

In (13.51) only the momentum of the scattered particle is resolved. The limit 
t — > oo is needed, so that for any initial condition on the plane t the particle has 
escaped from the scattering region. 

We rewrite (13.51) in a slightly more convenient form by splitting the generator as 

L = L + L! (3.6) 
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with 



Lxf(q,p) = p{q) / dp'2it8{oo(p) - uo{p')){2ir)-^d{p - p')f(q,p f ) . (3.7) 



Then 



e Lt = e L < 



ot + [ dse LoS L ie Loit - s) + f ds [ ds'e LoS ' Lxe^-^Lxe^-^ . (3.8) 
Jo Jo Jo 



Let 

v {p) = J dp'2it5{uj{p) -u(p'))(2ir)- 3/2 9(p-p') 
be the total cross section at momentum p and set 

fk(v,P) =S{p-k) exp - u(p) / dtp(q + Vu(p)t) 



(3.9) 



fk = S(p-k) exp - i/(p) / dtp(q + Vu(p)t) 



(3.10) 



Inserting ( 13. 8ft in ( 13. 5ft . the first term does not contribute, since k ^ k', and the 
second and third term have limits, since suppp is bounded. Then 

POO 

2n5(iu(k) - w(*0M*, k') = (/», Lx/+) + / dt</ fc " Lie L %/+> (3.11) 

./o 

with (-,-) denoting the scalar product in L 2 (M 3 x M. 3 ,dqdp). The time integral 
converges, since the probability to stay inside suppp decays exponentially. 



4 Summation over ladder diagrams 

The goal of this section is to establish that for u(k) = u(k') it holds that 



lim X 2 (27ryi lad (k',k;X)=a B (k,k > ). 

X — ^0 



(4.i; 



The prefactor A 2 balances the cross section of the scattering region which is of order 
A" 2 . 

We set 



and note that 



H±(h) =TiJ dk^y^eih - k 2 )G E ±(k 2 ) 
H + (k)+H_(k) = v{k) . 



(4.2) 
(4.3) 



The gate > 1 1 ■ corresponds to the integral operator 

ki k 2 k 3 

dk 2 G E+ (k 1 )G E+ (k 2 )X(2n)- 3 (V(k 1 - k 2 )V(k 2 - k 3 )) v 

= X(2n)- 3 G E+ (k 1 ) J dk20(ki-k2)G E+ (k 2 )px(h-k 3 ) 
= i\{27r)- 3 / 2 G E+ (k 1 )H + (k 1 )p x (ki - k 3 ) . (4.4) 

Therefore the effective medium propagator, k\ k, is represented by the kernel 



oo „ 

F_(Jfei, k; A) = 5(k 1 ~k) + J2 / dk 2 ...dk 

m=2 J 



J] {G^kM^r^H^Xpxik, - k j+1 )} , (4.5) 

3=1 



where in ( 14.51) we set k m+ i = k. We define 



Io{x; X) = I dkiexp[iX 1 (ki — k) ■ x]F-.(ki, k; X) 

dk ie ikl - x X 3 F_(k + Xh,k;X) . (4.6) 



Shifting the other internal momenta by k and rescaling by A yields 

oo „ m 

I (x; A) = f + dh... dk m e ikl - x Yl {XG E+ {k + A^)i(27r)- 3 / 2 

m=l j=l 

H+ik + Xk^pikj-kj^)}, (4.7) 

where in ( 14. 7p we set k m+ i = 0. H + is a smooth function. Thus H + {k + Xkj) may 
be replaced by H + {k). The Green's function G E+ has the integral representation 



XG E+ (k + Xk t ) = iX / dt^t^ik+xk^ (fc)-fe) 
Jo 

POO 

i i dte- it{yu{kykl - ie) . (4.8) 







Inserting in (14.7ft yields 
lim Iq(x\ A) 

OO p p 

1 + 53 ( - H + (k)) m dki... dk m (2ir)- 3m dx 1 . . . dx m p(x x ) . . . p(x m ) 

m=l 

/■OO /"OO r m 

x / dt\ . . . dt m exp — z 53 feVw(fc) • fcj + (A; 3 - — %+i) • Xj) + zfci 



i=i 



x /"OO /-c 

m=l - 70 - 70 



cxp 



xp(x - Va;(fc)ti) . . . p(x - Vu;(fc)(*i + . . . + t m )) 



-H+(k) / dtp(x + Vu(k)t) 



(4.9) 



With this input we compute the single collision diagram with effective medium 
propagators, 



k' 



k 



It is given by 



k' ' ' - k 



h{\) = J dk 1 dk 2 dg 1 dg 2 X(27ry 3 (V(h - k 2 )V( 9l - g 2 )*) v 

xF+(h, k'; X)F4k 2 , k; X)F+( 9l , k'; X)*F_(g 2} k; X)* 

= (27r)- 3 XX- 3 J dk 1 dk 2 dg 1 dg 2 (2n)- 3/2 J dxp{x)6{k' — k + Xfa — Xk 2 ) 

xe- ikl - x X 3 F„(k' + Xh, k'\ X) (e- i9l X X 3 F^(k' + Xg u k'; A))* 
xe ik2 - x X 3 F + {k + Xk 2 , k; A) (e i92X X 3 F + {k + Xg 2 , k; A))* . (4.10) 

Here F + is obtained from F_ by substituting px(—kj + for p\(kj — kj +1 ) in 

(14.5p . Inserting f]4.9f) one obtains 



A (A) 



(2tt)- 3 A~ 2 y c/xp(x)(27r)- 3/2 0(A; - k') 

- (H + {k) + H_(k)) / dtp(x + Vu;(A;)t) 



exp 



{H+(k') + H_(kf)) / dtp(x + Vuj(k')t) (l + 0(A)) . (4.11) 
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Using ( 14.31) yields indeed 

lim 2vr5(a;(A;) - oo(k'))X 2 (2n) 3 h(X) = (/ fc " LJ+) . (4.12) 

A — ^0 

As next item we consider the ladder diagram with two collisions, 

k' k 



— => — < 


< f-^» — 


— < — i 


( 1 — < — 



k' -^4 i-e— A; 



The external legs carry only the free propagator. Their gates can be summed to the 
effective medium propagator as in the first part of this section. The top line has 

m + 2 vertices and internal momenta ki, . . . , k m+ i, m — 0, 1, 2, Correspondingly, 

the bottom line has n + 2 vertices and internal momenta g 1 , . . . , g n+ i, n — 0, 1, 2, . . .. 
The sum over all gates reads 

oo oo 



J 2( A ) = J2J2 dkl '-- dk ™+ 1 d 9i--- dg n+ i0(k' - h)6(k m+1 - k) (4.13) 

m=0 n=0 ^ ^ 

x(27r) _6 Ap A (-/ci + gi)Xp x (-k m+1 + g n+1 )* 

rn 

J] {A(2tt)- 3 / 2 GW(M^4Mpa(% - 



X 

x J] {A(27r)- 3 / 2 G s _fe)(-0iJ_fe)p A fe - g j+1 y}G E+ ( k m+1 )G E „(g n+1 ] 

3=1 



We shift kj by fc m+ i and gj by g n+ i and then substitute k m+ i = u + ^v, g n+ i = u— ~v . 
Finally v is rescaled to Xv , kj to Afcj, and gj to Ag,. Then ( I4.13P becomes 

OO OO „ „ „ 

^2 (A) = / <i/ci . . . dk m / dgi . . . dg n / dudv9(k' — Xki — u — \Xv) 

m=0 n=0 J J J 

x9(-k + u + ^Xv){27T)- 6 X 2 X^ 3 p{-k 1 + gi - v)p{-vf 



rn 

X 



\\ {X(2n)- 3/2 G E+ (u + \Xv + Xk^iH+iu + \Xv + A%)p(% - fc i+1 )} 



n 

x J] {A(2tt)- 3 / 2 G^( M - \Xv + \ 9j )(-i)H-.(u - \Xv + A^)p(^ - 
xG^^ + ±Au)G E _(u - |Au) , (4.14) 
where in (I4.14p we set fc m+ i = , g n+ i = 0. 
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We insert the approximation (14.81) and switch to position space, 

OO OO p p p 

= /J / dk\ . . . dk m / dgi . . . dg n I dudv9{k' — u)9(u — k) 

m=0 n=0 J J J 

poo POO 

x( — H + {u)) ( — H_(u)) / dt\...dt m I ds\...ds n 

Jo Jo 

m 

x exp |^ — i'^2' t j(^ 1 ( LJ ( u ) — E — is) + Vuj(u) ■ kj + \Vuo(u) ■ v) 



+i^2sj(X 1 {u>{u) - E + ie) + Vw(«) • gj - \Vuj(u) ■ v) 

3=1 



x j dx x . . . Gte m (27r) 3m p{x x ) . . . p(x m ) j d Vl . . . dy n (2ir) 3n p( yi ) . . . p(y n ) 

r m 

x(27r)~ 6 A 2 A~ 3 (27r)~ 3 / dxdyp(x)p(y)exp -i^^Xj • (kj - k j+l ) 



3=1 



+ i y i ' (to ~ g i+i) ~ ix • +gi-v)-iyv 



3=1 



xG E+ (ti + \\v)G E -{u - \\v) (1 + 0(A)) 



(4.15) 



where in (14. 15f) we set k m +i = 0, g n +i = 0. We integrate over k±, . . . , k m} gi, . . . , g n 
and substitute 



(4.16) 



i=l 



Then 



OO OO „ 

7 2 (A) = J2J2 dudv{27i)- 3 6{k' - u)6(u -k)(- H + {u)) m ( - H_(u)) n 



m=0 n=0 

x / dT X ...dT m \ da x . . . da n (2iry 6 \ 2 \~ 3 

J0<Tl...<T m J 0<(71...<lT ri 

m n 

dxdyp(x)p(y) JJ p(x - Vu(u)Tj) JJ p(y - Vuj(u)(Tj) 

3=1 3=1 

x exp - ^(n))A _1 (r m - cx n ) - i\Vu(u) ■ v(r m + a n ) + i(x - y) ■ v ] 
xG E+ (u + \\v)G E -{u - \\v)(l + 0{\)) . (4.17) 

We still have to take the limit A — > for the two remaining Green's functions. 
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Collecting all the u, v dependence into the smooth function /, this leads to 

/3(A) = A / dudvf(u,v) exp [i[E - u{u))\~ l {T m - a n ) - \iVu{u) ■ v(r m + a n )] 



xGe+(u + \\v)Ge-{u — \Xv 



2' 

00 roo 



= A / dudvf(u,v) / dt ds 
J Jo Jo 

x exp \i{E - uj(u))\~ l (T m - a n ) - \iS7uj{u) ■ v(r m + a n ) 

xi(t -s)(E- u{u)) - e(t + s) - \iVu{u) ■ v\{t + s)] (l + 0(A)) 

/roo roo 
dudvf{u,v) I dt I ds exp [i(E — cu(u))(t — s) 

-liVu(u)-v(t + s)](l + 0(X)) . (4.18) 
We rotate s, t by 7r/4 and rescale t + s by A. This yields 

/roo 
dudvf(u,v) / dte- iVuj{uyvt 2ir5(E - lu(u)) . (4.19) 

</max(T m ,cr n ) 

Inserting (14.191) in (14.171) one extends the t- integration from to 00, while r m < t, 
a„ < t. Then 



f 00 

J 2 (A) = (2vr)~ 3 A- 2 J dt J dudv(2Try 3 e\k' - u)6{u - k)(2ir)- :i 

dxdy p(x) p(y)2n5 (E — u>(u)) exp [ — iVuj{u) ■ vt + i(x — y) ■ v\ 



x exp — H + (u) / dsp{x — Vuj{u)s) — H_{u) I dsp{y — Vu{u)s) 
L Jo Jo 

x(l + C(A)). (4.20) 
Let us set = S(p — k). Then (I4.20p amounts to 

roo 

lim 2ir5(uj(k) - u(k'))\ 2 (2Tr) 3 I 2 (\) = / dt(f kl L^Ljj!,) . (4.21) 

Summing over the gates in the external legs amounts to replacing by f£ 
and fk by fj7. The ladder diagrams with increasing number of rungs generate the 
time-dependent perturbation theory for e Lt . Thus we have verified the limit (14. ip . 

5 Backscattering and maximally crossed diagrams 

The first subleading contribution to (<J\(k, k'))y comes from the maximally crossed 
diagrams. They are responsible for a narrow peak of width A centered at k! = — k 
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as a correction to ob(A;, k'). Recall that (fk(q) = (27r) 3 ' 2 e tk ' g . We define 

2tt,5(^(A;) - ^(A;'))a B (fc, fc'; k) = (2vr) 3 / dt(/ fc "^, L ie Lt L lV J+) (5.1) 

</ o 

and 

o-back,fc(«) = <r B (k t -k; k) . (5.2) 

Note that the first and last scattering event carries now an extra phase factor. If 
one chooses u(k) = u(—k + Xk), then the backscattering rate is given by 

lim(27r) 3 A 2 J max (A;, —k + Xk; X) = a hBck}k (K) . (5.3) 

A — >0 

( I5.3p is our main result. 

To verify our claim does not require a long computation. Let A max be a particular 
maximally crossed diagram. If its internal momenta at the lower line are g%, . . . , g m , 
then we substitute —g m -j for g m . Then A max (k', k) = A\ ad (k', k) where A\^{k\ k) is 
the corresponding ladder diagram with external momenta k! and k for the upper line 
and —k and —k' for the lower line. In particular A max (—k, k) = Ai a a{— k, k), which 
implies ( 11. 3p . I.e. for precise backscattering at small A the scattering rate is twice 
the one predicted by the Boltzmann equation, upon omitting the single scattering 
contribution. This result is independent of the shape function p. 

For the fine structure we look Xk away from backscattering and have to sum over 
all ladder diagrams, as explained in Section 4, with the correspondingly modified 
external momenta. 

Let k\,k2, resp. gi,g% be the internal momenta of the last rung of the ladder 
diagram. Then the diagram beyond the last rung equals 



I r (k, -k + XK; A) = J dk 2 dg 2 X(V(k 1 -k 2 )V(g 1 -g 2 )*) v F_(k 2 , k-Xn; A)F_(^ 2 , k; A)* . 

(5.4) 

We shift k 2 to k — Xk + k 2 and g 2 to k + g 2 and rescale, 

I r (k, —k + Xk; A) = AA~ 3 J dk 2 dg 2 Q{k\ — k — Xk 2 + A/t)p(A _1 (A; 1 — g{) — k 2 + g 2 + k) 
xA 3 F_(/c - Xk + Xk 2 , k - Xk; A)A 3 F„(A; + Xg 2 , k; A)* . (5.5) 

Thus, compared to the ladder diagram, in position space the diagram picks up the 
phase factor e~ iK ' x . Repeating the argument for the leftmost rung one concludes the 
validity of (I5.3p . 



6 Backscattering in the diffusive approximation 

Let us consider the conventional slab geometry for which p{r) = 1 if r 3 > and 
p(r) = otherwise, r = (ri,r 2 ,rs), and, for the sake of illustration, choose k = 
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(0, 0, fc 3 ), &3 > 0. The scattering rate is proportional to the cross sectional area, 
through which we divide by imposing that the first scattering is at r = (0,0, r 3 ) 
with r 3 > 0. Even with these simplifications the inverse of L, as needed in ( 15.11) . 
is not so easily computed. Therefore we approximate the motion between the first 
and last scattering through a Brownian motion, i.e. L through D(u)A r , where 
D(u) = D(u(k)) is the diffusion coefficient obtained from the transport equation 
(11. ip . According to this equation the particle simply leaves the scattering region 
upon hitting its boundary. Thus A r is taken with Dirichlet boundary conditions at 
{r 3 = 0}. Setting G D = (A r ) _1 one has 

G D (r;r') = (4rt\r -r'])- 1 - (4rt\r -r'])- 1 , r 3 ,r 3 >0, (6.1) 

where r — (ri,T2, — r 3 ). The first and last scattering are approximately a mean free 
path away from the boundary {r 3 = 0}. Therefore the diffusion approximation is 
taken between the points r = (0,0, £*) and r' = (r\\,£*) with the mean free path 

t* = \V k Lu(k)\/v(k). 

The backscattering in the diffusive approximation is thus given by 

Cback,D(^ + Aft) = (K^ 

) 2 /D(uj)) I dr||G D (0,0,r;r||,r)(l + cos(«-r||)) (6.2) 

with u(u) = u(u(k)) and k = (ki,k 2 ,0). In the round bracket to the right the 
"1" accounts for the incoherent scattering while the "cos" expresses the coherent 
backscattering. Working out the integral one obtains 

(T hacKD (k + A«) = (v(u) 2 /D(u))t(l + (2t\n\)- 1 (l - exp[-2^*|«|])) . (6.3) 

We refer to [6] for a confirmation of (16. 3p in a light scattering experiment. 



7 Light scattering 

As an illustration we consider the wave equation with a random index of refraction. 
Other wave equations can be handled in a similar fashion. The wave field <fi : R 3 — ► R 
is governed by 

— (j) = c(x) 2 A ( p (7.1) 

and we set 

c( x ) = 1 + y/XV(x) (7.2) 

with V from (12 .4p . The Gaussian statistics looks unphysical because c(x) will have 
regions where it is negative. However for a, say, bounded random potential only the 
Gaussian part persists for small A [3| [TT]. 

The general strategy is to rewrite (17.11) in the form of a Schrodinger equation 
and then to apply the results from before. Such a procedure is not unique and we 
adopt the one employed in [11], see also [12]. 
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The dispersion relation for (17. II) at A = is 

w(*) = |*|. (7.3) 

Again, more general dispersions can be handled by our method. If we consider uo as 
multiplication in Fourier space, then the corresponding operator in position space 
is denoted by Q. We introduce the two-component "wave function" 

V> = (^~), ( 7 - 4 ) 

where 

tp ± {x) = -=(n<f)(x)±i(l + y/XVix))' 1 ^)) . (7.5) 



Then (17. ip is rewritten as 

d A/> + \ ftt \ fip + \ rr ( VVL + VLV VVL-VlV \ (i). 



i- 



^ A T/n , nt/ T/n ni/ ./. ■ ( 7 - 6 ) 



dt\?P-J v° -ttJvl>-J \-vn + nv -vvl-vlv)\^ 

Note that ||^+|| 2 = ||V>~|| 2 is the energy of the wave field, hence conserved in time. 
We regard ( 17.6P as an evolution equation in L 2 (IR 3 ) ©L 2 (M 3 ). The physical solutions 
are then the subspace defined by (17.51) . The right-hand side of (17.6p defines a self- 
adjoint operator, thus the solution conserves the norm. 

The only still missing input are the plane wave solution in case A = 0. They 
are of the form <f>k{x) = {2n)^' i / 2 uj{k)^ 1 exp[i(x • k — u(k)t)], where the prefactor is 
chosen such that the 0fc's are orthogonal in the energy norm. Thus 

tfj + {x) = {2-K)- z ' 2 e lk - x , iff~(x) = . (7.7) 

We have 

(VttijIXk) = (2tt)-^ 2 J dhV{k - ki)u(ki)ip(ki) , 

(QV^jJ{k) = {2Tx)-' i ' 2 uj(k) [ dk x V{k - k^kx) . (7.8) 



In the perturbative expansion, compared to the Schrodinger case, a vertex carries 
an extra factor of u. This amounts to an extra factor of uj 2 in the collision rate. 
Thus the transport equation for wave propagation has the backward generator 

L w f(q, P ) = Vw(p) • V q f(q,p)+p(q)uj(p) 2 J dp'2^5{uj{p) - u(p')) 

x(2n)- 3 / 2 e(p-p')(f(q,p')-f(q,p)) , (7.9) 

compare with (13. 2p . 

For the coherent backscattering of light one simply replaces in (15. ip L by L w . In 
other words the collision rate v(p) is replaced by oj(p) 2 v(p) and the propagation is 
Vu(k) = k/\k\ rather than k. In the diffusion approximation such fine details are 
no longer visible. Provided i* is taken as the mean free path for light, the shape of 
the coherent backscattering peak remains unaltered. 
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